INFORMATION LENGTH AND LOCALIZATION IN ONE DIMENSION 

Imre Varga and Janos Pipek 

Quantum Theory Group, Institute of Physics, 
Technical University of Budapest, H-1521 Budapest, Hungary 

*rf ' Abstract - The scaling properties of the wave functions in finite samples 

Q\ . of the one dimensional Anderson model are analyzed. The states have been 






PACS. 71.50 - Localized single-particle electronic states 
PACS. 71. 55 J - Localization in disordered structures 
PACS. 72.15R - Quantum localization 



characterized using a new form of the information or entropic length, and 
compared with analytical results obtained by assuming an exponential envelope 
function. A perfect agreement is obtained already for systems of 10 3 -10 4 sites 
over a very wide range of disorder parameter 10 -4 < W < 10 4 . Implications 
for higher dimensions are also presented. 
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The numerical detection of exponential localization in finite, random sys- 
tems is not a trivial task, especially in the weak disorder limit when the local- 
ization length of the eigenstates is expected to be larger than the system size. 

•i-H . 

^ This may be a problem even in one dimension (ID), where rigorous results [1] 

affirm complete exponential localization for any strength of disorder. 

The model under consideration can be described by a tight-binding Schro- 
dinger equation in the nearest neighbor approximation as 

U n +1 + Un-1 + V n U n = Eu n , (1) 

where V n are independent random variables with uniform distribution over 
the [— W/2 . . . W/2] interval, E is the eigenenergy and u n is the amplitude 
of the eigenfunction on site n. The wave functions are expected to behave 
asymptotically up to oscillations as 

u n ~ exp(-7n), (2) 
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where 7 = £ 1 is the inverse localization length or Lyapunov exponent that 
may be numerically obtained as [1] 
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The exponential localization has been corroborated by the one-parameter 
scaling theory [2]. This was built essentially on the basis of the Thouless 
number [3] that is related to the dimensionless conductance [4]. Inspired by 
that scaling law in a recent paper Casati et al. [5] have introduced the concept 
of information localization length [6] and investigated numerically the possible 
scaling properties of the eigenstates themselves. Their study has also been 
motivated by previous results obtaining a scaling law for an analogous model 
in quantum chaos, the kicked rotator [7], as well as for band random matrices 
[8] using the same concept. 

The main idea of Ref. [5] is to take a suitable ensemble of eigenstates and 
to calculate the information length as 

p c (E,N,W) = eM~S-Sref), (4) 

where N stands for the system size, W describes the strength of disorder. Index 
C is used to label (3 in order to refer to the definition of Casati et al. S is the 
averaged information entropy, 

N 

S = -J2 u l lnu2 m ( 5 ) 

n=l 

of normalized eigenstates in a window around energy E for different realizations 
of the random potential. The S re f in Eq. (1) stands for the entropy of a 
reference state 

u n ~ sm(<pn) where cos<p — E/2. (6) 

This wave function is the exact solution of Eq. (1) in the absence of disorder 
(W = 0) with uq = and u\ — 1. A straightforward calculation yields the 
asymptotic form S re f(N) — ► ln(2Af) — 1 as iV — ► 00. We have to indicate that 



the particular choice of S re f in Eq. (4) involves a delicate problem, that we 
wish to discuss below. 

The principal aim of definition (4) is that (3c should give the portion of 
the sites significantly populated by the eigenstates compared to that of the 
reference state (6). It is clear that with the increase of the system size iV we 
expect exp(S') oc N in the case of extended states and exp(S') — > const, for 
localized states. Note, that since for large N, exp(—S re f) ~ c/N, its role is 
normalization. 

Casati et al. have numerically established the scaling law [5] 

with C ~ 1. ^oo has been calculated as ^ = 1/7jv (cf. Eq. (3)) for strong 
disorder (1/7jv < N). For weak disorder (1/7jv > N) instead of numerically 
calculating the Lyapunov exponent the authors of Ref. [5] used £oo given by 
the perturbative calculation [9]. They have found, however, no theoretical 
explanation for (7). 

We would like to point out that the approach of Casati et al. [5] resides 
on the supposition that in Eq. (7) (3c < 1 i.e. S < S re f- There is, however, 
no rigorous proof for these relations. This is particularly crucial in the weak 
localization limit W — ► 0, (3c ~ 1- Indeed, in actual numerical calculations for 
small ensembles we have found occasionally (3c > 1, as well. We could prove, 
however, that S < S re f (N) for sufficiently large ensembles. The proof is based 
on the idea that any small perturbation {Q' n = Q n (l + S n )}, 5 n = of an 
arbitrary probability distribution {Q n } leads to S < S, where S is an average 
over an infinite number of realizations of S n . As a byproduct, it is also evident 
that for the calculation of (3c, S re f = S re f(N) has to be applied instead of 
using the asymptotic form ln(2iV) — 1. 

Furthermore, as we will show it later, for finite systems the functional form 
of the scaling law (7) does not seem to be valid in the weak localization limit 
(W — > 0). This is mainly due to the improper application of the perturbative 
treatment. 



In this Letter we propose a new form of (3 for resolving the above men- 
tioned difficulties and appropriately handling the W — »■ limit. We will demon- 
strate that our (3 function is, apart from showing special scaling properties, 
capable to prove the exponential localization in one dimension. 

On the basis of our recently introduced classification scheme [10], the 
information entropy (5) of any general, normalized, nonnegative lattice distri- 
bution can be split as a sum of two terms 

S = S str + lnD, (8) 

where D is the derealization measure [11] or participation number [12] 

and S str is the structural entropy of the distribution. Parameter D is widely 
used in the literature giving the number of sites the eigenstate extends to. 
Therefore it is bounded as 1 < D < N. Using D we may introduce a normalized 
quantity q, the spatial filling factor or participation ratio as 

q = D/N, for which < q < 1. (10) 

The structural entropy in equation (8) has been shown [10] to be nonnegative 
with bounds 

0<S str <-lnq. (11) 

Using the quantities discussed above we now propose an alternative form 
of the normalized information length of an eigenstate as 

/3 = ^expS, (12) 

which using expressions (8) (9) and (10) becomes 

(3 = qexp(S s tr)- (13) 

Due to the well-known properties < S < IniV the following restrictions are 
imposed on (3 

< (3 < 1. (14) 



These bounds are valid for each state separately, whereas fie < 1 can only be 
guaranteed after an appropriate averaging process. Obviously, an eigenstate 
expanding uniformly over the whole system will have v? n = 1/N so that (3—1. 
In the other extreme, for localized states D ~ 1 and S s t r ~ 0, in a finite system 
one obtains (3 ~ 1/N. 

Returning now to the Anderson model (1) with non-zero disorder W ^ 0, 
one expects according to (2) and (6) the charge distribution of the solution of 
the form 

\ u n\ 2 ~ fiin) sin 2 (an), (15) 

where / is a slowly varying envelope function due to the presence of the per- 
turbing random potential. Obviously, in our case function / is expected to take 
an exponential form f(p) = exp(— p). The value of a is roughly ip defined in 
(6). Further specification of a is needless for our purposes, the only restriction 
we impose is a ^> 1/N that is always fulfilled except very close to the band- 
edges. As we have shown [10,13] for multiplicative superstructures of the form 
(15) we get 

In q = In q* + In q (16) 

and 

£>str = J s tr + '-'sir' (1 ' ) 

where the upper index / stands for the values obtained for the charge distri- 
bution f{^n) alone and upper index indicates the ones for sin 2 (an). It is 
possible to show that (independently of a) 

q° = 2/3 and S° str = In 3 - 1. (18) 

Using (16) and (17) as well as definition (13) we get 

[3 = [3 f [3o, (19) 

where (3q = q° exp(S® tr ) ~ 0.7357. In the limit of vanishing disorder W — > 
one expects 7 — ► 0, f(p) — ► 1. For such distribution using (5), (8) and (9) 
qf -> 1 and S f str - therefore f3 f - 1. For strong disorder on the other hand 



7 >> 1 which yields D ps 1, therefore in finite systems qf ~ 1/iV and S{, tr ~ 
resulting /3f « 1/iV. 

The role of /?o in (19) is similar to the factor exp(SVe/) in (4) introduced 
by Casati et al. [5] , however, (3q is a constant independent of the system size 
N, and its derivation is based on the separation of the wave function to an 
envelope and a strongly oscillating part (15). 

Before performing the numerical simulation we still have to give the ex- 
plicit qf and SJ, tr values as a function of 7. In Ref. [10] the general form of qf(z) 
and S str (z) functions with z = 7./V = JV/£ is given for arbitrary dimensionality, 
that have been calculated applying a continuous lattice approximation, i.e. the 
relevant scale (7 _1 ) was assumed to extend over many lattice spacings. It has 
been shown in [10] that in one dimension 

*«) = §f , (20a) 

and 

s '^ = W)Hm)' <m) 

where functions F(z), G(z), and H(z) are defined as 

F(z)= [ f(p)dp, (21a) 

./o 

G(z)= f f(p)dp, (21b) 

Jo 

H(z) = - [ Z f(p)Hf(p)]dp. (21c) 

Jo 

Inserting the general functions given in expressions (21) into (20) one obtains 
for f3f in the continuous limit 



/*/(*) 



qf(z) eMSLi*)) = z- 1 F(z) exp (?^j . (22) 



It is straightforward to calculate the (3f(z) function for any envelope shape 
f{p). For exponential decay, f(p) = exp(— p), we get 



exp z — 1 

z exp z \ exp z — 1 



/W^) = ,1 . ex P ( l " — .. ; ) • ( 23 ) 
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Expressions (22) and (23) are the principal results of this Letter showing the 
scaling property of j3f provided that a reasonable definition for the f(p) decay 
function (15) exists. 

Let us turn now to the asymptotic properties of (3/(z). In the case of 
strong localization z — ► oo (N ^> £), since both F(oo) and H(oo) are finite for 
most of the practical cases, from the general expression (22) one gets 

/3 / (^)~^ 1 ~|, (24) 

as expected. In the other limit of derealization as z — ► (e.g. £ — ► oo keeping 
N fixed) we found that the asymptotic form of f3f (z) is governed by the short 
range properties of the form function f(p), i.e. it depends on the derivative of 
f(p) at the origin p = 0. Namely, if /'(0) ^ then 

(3(z) « 1 - ^z\ (25) 

while for a Gaussian form function e.g. where /'(0) = the first nonvanishing 
term is of the order of z A . 

Instead of the flf(z), after Ref. [5], we define 

in order to emphasize both the localized and delocalized limits. In our numer- 
ical simulation we have compared Eq. (26) for exponential form function (22) 
with the calculated 

-T^ (2?) 

values, where (3/ is defined here as 

(3 f = —qexpS str - (28) 

Po 

As < f3f < 1 is true already for individual wave functions it is not neces- 
sary to perform averaging over states taken from an energy window, however, 
at a certain energy we calculate the statistical means q and S s t r over many 
realizations of the random potential. 



The eigenvectors in our simulation were obtained by the iteration of the 
recurrence relation of Eq. (1) with initial conditions uq = and u\ — 1. In all 
of the presented results the length of the system was N = 10 4 and the number 
of samples used for averaging was M = 10 3 . The energy was fixed to E = 0.1. 
The localization length was obtained as 

Z(E)=-?(E)-\ (29) 

were 'j(E) was calculated according to Eq. (3). In a finite lattice two relevant 
length scales characterize the system the chain length N and the lattice spacing 
(a = 1). The relation of £ with respect to these length scales is essential in 
such type of calculations. 

In Figure 1 we have plotted the localization length as a function of the 
strength of disorder ranging from W = 10 -4 up to W = 10 4 . Our numerical 
calculations confirm the theoretically expected behavior £ _1 ~ InW for large 
disorder. In the case of vanishing disorder W — > perturbation theory pre- 
dicts £ _1 ~ W 2 in the thermodynamic limit N — *■ oo [9]. This consideration, 
however, fails for finite N as one can see in the low W part of Figure 1, where 
we have numerically found £ _1 ~ W. The figure clearly shows that this finite 
size behavior becomes relevant for such disorder values where the localization 
length is comparable to or greater than the system size. It seems that the 
behavior of the wave function on intermediate length scales is governed by a 
different characteristic length £$. As is shown in Figure 2 on this latter scale 
exponential localization can be detected, as well. 

In Figure 2 we have compared ftc of Casati et al. and our numerical 
f3 f and theoretical /9 exp as a function of £/N. In each of these cases we have 
plotted y = (3/(1 — (3) vs £/N in a log-log plot. The results of the simulation 
follow the curve for exponential localization. The deviation for ln(£/iV) < —9 
is present because the value of the localization length becomes comparable 
to the lattice constant as W — ► oo i.e. we obtain Iny — >■ — IniV. Both our 
analytical and numerical results confirm the expected high disorder behavior 
of lny m In/? ~ ln£ (see e.g. (24)), as well as the low disorder behavior of 
Iny pa ln(l — (3) ~ 2 ln£ (see e.g. (25)). For strong disorder we observe perfect 
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agreement with the scaling law set up by Casati et al. [5] in Eq. (7). Our 
analytical results (23) confirm this scaling law 

\n(3 exp (z) ->-lnz + l, (30) 

in the limit of strong localization z = iV/£ — ► oo. 

For weak disorder, however, Figure 2 shows a considerable disagreement 
between our results and that of Ref. [5] . This is easy to understand considering 
that for this regime Casati et al. have used the predictions of the perturbative 
calculation valid in the thermodynamic limit. As we pointed out earlier this 
approach needs a careful analysis. They compare quantities y and £oc where 
y is calculated from wave functions characterized by the intermediate length 
scale £j. On the other hand calculating the £ according to Eqs. (3) and (29) we 
obtain an almost perfect agreement between the numerical simulation and our 
analytical expressions for exponential form function. This shows that, apart 
from the exponential long range behavior for ^ < N, in the intermediate range 
(£oo > N) the same kind of decay was found with a different scale constant £$, 
as well. 

Just to have a feeling how well the charge distribution of the form of Eq. 
(15) describes the average properties of the wave functions in the Anderson 
model, in Figure 3 we have plotted S s t r as a function of q. We have compared 
the results of the simulation with analytical results obtained assuming the 
form of Eq. (15). A satisfactory agreement can be established between the 
numerical and analytical results especially for low and high disorder. The 
charge distribution is clearly not of pure exponential form, but a plane wave 
modulated by an exponential envelope. 

We would like to note that especially for weak disorder, the exponential 
localization, apart form fluctuations and oscillations, is still strictly true and 
visible using our construction of /3-function at least up to localization lengths 
several times larger than the size of the system. We believe that after a proper 
definition of the localization length, a similar procedure could clearly show the 
expected exponential localization in two dimensions, as well. Results along 
these lines are to be published in a subsequent paper. 
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Figure Captions 

Figure 1. The log- log plot of the localization length versus strength 
of disorder using the numerical simulation. The characteris- 
tic relations are also denoted. 

Figure 2. Scaling of the entropic length versus the localization 
length using \ny = ln(/3/(l — (3)) versus ln(£/iV). Solid 
symbols represent the results of our numerical simulation. 
Dashed line stands for the scaling law found numerically by 
Casati et al. [5] . Note the change in the slope of the contin- 
uous curve at around £ ~ N. The continuous curve is our 
analytical result (see Eqs. (23) and (26)). 

Figure 3. The structural entropy versus filling factor in a semi-log 
plot. The filled symbols represent our simulation, while the 
solid line stands for the relation assuming a charge distribu- 
tion of the form of Eq. (15). 
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